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ABSTRACT 

We suggest a revised distance to the supernova remnant (SNR) G109. 1-1.0 (CTB 
109) and its associated anomalous X-ray pulsar (AXP) IE 2259+586 by analyzing 
21cm Hl-line and 12 CO-line spectra of CTB 109, HII region Sh 152, and the adjacent 
molecular cloud complex. CTB 109 has been established to be interacting with a 
large molecular cloud (recession velocity at v = —55 km s _1 ). The highest radial 
velocities of absorption features towards CTB 109 (—56 km s _1 ) and Sh 152 (—65 
km s _1 ) are larger than the recombination line velocity (—50 km s -1 ) of Sh 152 
demonstrating the velocity reversal within the Perseus arm. The molecular cloud has 
cold HI column density large enough to produce HI self-absorption (HISA) and HI 
narrow self-absorption (HINSA) if it was at the near side of the velocity reversal. 
Absence of both HISA and HINSA indicates that the cloud is at the far side of the 
velocity reversal within the Perseus Arm, so we obtain a distance for CTB 109 of 4±0.8 
kpc. The new distance still leads to a normal explosion energy for CTB 109/ AXP IE 
2259+586. 

Key words: supernova remnants: individual (CTB 109)-HII regions: individual (Sh 
152), pulsarsdndividual (AXP IE 2259+E 



1 INTRODUCTION 

The Supernova Remnant (SNR) CTB 109 hosts an anoma- 
' lous X-ray pulsar (AXP) IE 2259+586 with a superstrong 
magnetic field, i.e. a magnetar. IE 2259+586 e mitted a Soft 
Gam ma-ray Repeater-like outburst in 2002 (|Kaspi et al.l 
2003), so CTB 109 is an especially interesting target for 
further study. CTB 109 is believed to strongly interact with 
the large molecular cloud (£=108.77°, 6=-0.95°) adjacent to 
its west si de, based on its semi-circular shape in b oth X-rays 
and radio (|Tatematsu et al.lll99d : ICoe et al.ll 19891 ). CTB 109 
has an extended X-r ay bright region in the interior, which 
is thermal in origin (|Rho fc Petrel 1 19971 ) and has no mor- 
phological correlation with AXP IE 2259+586 jPatel et all 
12001]) . The lobe likely originates from an interaction of the 
SNR shock with a cloud at —55 km s -1 (jSasaki et al.ll2004 

l2ooeh . 

Distance measurement to a SNR/AXP sys- 
tem is extremely important but usually quite dif- 
ficult. Many efforts have previously been aimed 
at di stance dete r minati on of SNR/pulsar sys- 
tems l|Camilo et al l 120061 : iGotthelf fc Halpernl 120091 ; 
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Leahy fc Tianl 2008; McClurc-Griffith & Gacnslerl 120051; 
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2006 

(|2002l . K02 hereafter) analyzed HI ab- 
emission spectra of CTB 109 and near 



compact radio sources, and estimated a distance of 3.0 kpc 
to CTB 109. 

On the other hand, iDurant fc van Kerkwiik! (|2006l . 
DvK06 hereafter) used the red clump stars method to mea- 
sure distances to AXPs. This method allows deriving the 
function of reddening versus distance in any given line of 
sight, based on field stars over a relatively small area. DvK06 
excluded a distance of 3 kpc for AXP IE 2259+586 based on 
its high extinction, instead preferring a much larger value of 
6 kpc. 

The conflict between these two distances motivates us 
to re-analyze distance evidence for CTB 109 using 1420 
MHz continuum, 21 cm HI- line and 12 CO data. We pro- 
vide a revised distance to CTB 109 by employing improved 
HI spectral analysis methods, and by analyzing possible HI 
Narrow Self-Absorption (HINSA) and HI Self-Absorption 
(HISA) from the adjacent molecular cloud complex. Based 
on similar considerations, the methods have recently been 
used to measure kinematic distances to Galactic molecu- 
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Figure 1. Left: 1420 MHz continuum image of CTB 109 and the near HII region Sh 152. Right: CO emission at —51 km s 1 with contours 
(8, 15, 21, 29, 100 K) of the continuum emission. The left panel shows boxes (both green and red) used for extraction of HI spectra 
(see text) and centered at [I = 108.80° ,b = -0.96°] for box 1, [108.76° , -0.98°] for box 2, [109.43° , -1.36°] for box 3, [108.76°, -0.96°], 
[108.84°, -1.02°] for box 4, [108.77° , -0.95°] for Sh 152 and [108.93°, -1.22°] for CTB 109, respectively 



lar clouds ijRoman-Duval et alJliool iKrco et alJl2008h . We 
compare our new results to the previous conflicting distance 
estimates of K02 and DvK06. 



2 RESULTS AND ANALYSIS 

We show the 1420 MHz continuum image (left) and the CO 
image at v = —51 km s~ (right) in Figure 1, both centered 
on CTB 109. As suggested by previous studies, there is a 
good correlation between CTB 109 and CO features in the 
velocity range -50 to -57 km s" 1 (see Fig. 4 of K02), but 
no respective convincing HI features in the velocity range 
associated with CTB 109. 

2.1 HI absorption 

We construct HI emission and absorption spectra for CTB 
109 and for Sh 152. Since CTB 109 is extended rather than 
a point source, the standard formula for HI absorption spec- 
tra does not apply. The continuum emission from CTB 109 
extends into our background region because the background 
region is chosen to be near the continuum peak. This min- 
imizes the difference in the HI distribution along the two 
lines o f sight (source and background, see iTian fc Leahvl 
|2008bl ). The resulting difference in brightness temperature 
is given by: AT = T off (v) - T on (v) = (T s a -TQ(l-e- r ^). 
Here T on (v) and T ff(v) are the average HI brightness tem- 
peratures of spectra from a selected area on a strong con- 
tinuum emission region of the source and an adjacent back- 
ground region which contains weaker continuum emission. 
Tg and T§ g are the average continuum brightness tempera- 
tures for the same regions respectively. t c (v) is the HI optical 
depth from the continuum source to the observer. 

The regions used for source and background spectra are 
shown by the solid-line and dashed-line boxes (red; the back- 
ground area excludes the source area) in the left panel of Fig. 
1. Fig. 2 gives the HI and 12 CO spectra. The spectrum from 
the south filament of CTB 109 (top in Fig. 2) shows a high- 
est HI absorption velocity of —56 km s _1 . The HI absorption 
spectrum from Sh 152 (botton) is the same shape as that 



given by K02, and shows a highest absorption velocity of 
—65 km s _1 . 

Sh 152 sits at the peak of a large high brightness- 
temperature CO molecular cloud at velocities of —48 to— 56 
km s _1 (see Fig. 4 of K02, and bottom in Fig. 2). There 
is clear HI absorption against Sh 152 in the same velocity 
range. This means cold HI gas surrounding the CO cloud is 
responsible for the HI absorption and is located in front of 
Sh 152. 

The CO spectrum of CTB 109 shows narrow CO emis- 
sion at —48 to —54 km s _1 in the direction of the southern 
filament. The fact that the highest radial velocities of HI 
absorption features towards both CTB 109 (—56 km s _1 ) 
and Sh 152 (—65 km s" 1 ) a re larger than th e recombination 
line velocity (—50 km s -1 . iLockmanl |l989) of Sh 152 con- 
firms the existence of a velocity reversal within the Perseus 
arm. This reversal explains why the HI gas at higher radial 
velocity is in front of Sh 152 which is at lower radial velocity. 

Based on Fig. 2, the HI column densities integrat- 
ing over all velocity in front of CTB 109 and Sh 152 are 
N^ ai ~ 4 and 6 x!0 19 T s cm" 2 respective ly (using N H / = 
1.82xl0 18 rA«r.. lDickev fc Lockmanl (|l990l )). For CTB 109, 
the column density of HI from —48 to —54 km s _1 is 
~3xl0 19 T s cm -2 (3xl0 21 cm -2 if T s ~ 100K). 



2.2 HI self-absorption and narrow self-absorption 

HI self-absorption (HISA) is generally produced by tem- 
perature fluctuations in the cold neutral medium (CNM). 
HISA features in the ne ar Perseus arm have been widely 
revealed in the CGPS (|Gibson et alj 120051 ) due to the 
above-mentio ned velo c ity re versal caused by spiral den- 
sity waves (|Roberts! 1 19721 ). A particular kind of HI 
self-absor ption, referred to as HI narrow self-absorption 
(HINSA, |Li fc Goldsmith! (|2003l )) is produced by cold HI 
mixed with molecular clouds and cooled by its molecular 
environment. HINSA is shown to share the spatial and kine- 
matic stru cture of cold mole cular clouds traced by CO emis- 
sion lines ( Krfio et al.ll2008T ). HINSA can only be detected 
in specific molecular cloud regions which have narrower line 
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Figure 2. HI and 12 CO spectra of CTB 109 and Sh 152 



width (1 to 2 km s _1 ) than background HI emission from 
CNM components. HINSA bears special relevance to this 
work because the originating location of the cold HI can 
be constrained by studying CO clouds at the same velocity 
given their co-location. 

The CO image and spectrum in the direction of Sh 152 
show a large molecular cloud complex in the velocity range of 
—48 to —56 km s _1 . Part of the CO cloud interacting with 
CTB 109 has been believed to be a major factor causing 
CTB 109 to have an asymmetric semi-circular shape struc- 
ture, i.e. a relatively bright eastern part, well-defined in radio 
and X-ray images, and no emission in the west. Interaction 
of the SNR shock with the part of this molecular cloud com - 
plex at —55 km s _1 has been observed l|Sasaki et al.l l2006). 
HINSA could be produced if these CO clouds are in the front 
of warm HI background. 

We searched the HI channel maps and HI emission spec- 
tra over the large region covered by the CO cloud complex, 
and found no HISA whose HI spectra show dips of 10 K 
with the FWHM Av of the absorption line 5 km s _1 . This 
is consistent with the non-detection of HISA features in the 
region in the CGPS HISA census (|Gibson et al.ll2005h . As 
an example, Fig. 3 shows the HI emission spectra and re- 
spective CO spectra extracted from boxes 1 to 4 shown in 



Figure 3. HI and CO emission spectra from boxes 1-4 



Fig. 1 (green). The top two rows of Fig. 3 clearly show no 
HISA features from the CO cloud at — 51±5 km s -1 against 
the warm HI background (Tjh~ 100K in the range of —40 
to —60 km s _1 ). In contrast, a HISA feature from box 3 
nearby CTB 109 is seen (the third row) . We found one area 
of the CO cloud complex showing narrow CO lines and still 
bright enough CO emission (at Z=108.85° and 6=-1.02° that 
is labelled "4" in Fig. 1) to satisfy the conditions needed to 
produce HINSA. The CO emission spectrum from box 4 
shown in the bottom of Fig. 3 has a narrow line width of 
2 km s _1 . The respective HI emission spectrum shows no 
HINSA features. 
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3 DISCUSSION AND CONCLUSION 

3.1 The velocity /distance relation within the 
Perseus arm 

21 cm HI observations can provide distance estimates to 
Galactic objects using the radial velocity-distance relation 
based on the flat circular rotation curves model. However, 
this rotation model may lead to serious overestimation of an 
objects' distance, especially at the Perseus arm where the 
velocity field i s strongly influenced by the spiral arm shock 
l|Robertsll 19721 ) . The velocity reversal inside the Perseus arm 
is caused by a spiral shock, while in other parts of the outer 
Galaxy the radial velocity decreases monotonically with dis- 
tance. This velocity reversal creates a distance ambiguity for 
objects with radial velocities of — 50 km s -1 to —60 km s" 1 
in the line-of-sight to CTB 109 l|Tian et all 120071 , TLF07 
hereafter) . 

W e use the two-arm HI model of lFoster fc MacWilliamsl 
J2006h (FM06 here after), which is an updated version of 
the iRobert si j 19721) model. FM06 fit an empirical model 
of Galactic structure and density-wave motions to observa- 
tions. A comparison between photometric and the HI model 
distances to 22 objects (19 HII regions and 3 SNRs), shows 
individual point errors of 0.3 to 0.8 kpc (Fig. 8 of FM06). 
The 20% error quoted in FM06 is dominated by systematic 
errors. The velocity of -55km s" 1 for CTB 109 gives 3.2 kpc 
in FM06 model if it is at the near side of the velocity reversal 
(Fig. 6 of TLF07). Since our new result is that CTB 109 is 
beyond the velocity reversal, the FM06 model gives 4.0 kpc 
distance. The difference in distance is 0.8 kpc. The error in 
this value is dominated by the systematic error of 0.16 kpc. 

Because of the velocity reversal in the Perseus arm, 
HISA and HINSA should be produced if the cold CO clouds 
associated with CTB 109 are in front of the warm HI back- 
ground, i.e. at distance of 3.2 kpc. Absence of any HISA and 
HINSA features indicate that the cold CO cloud is at the 
far distance of 4 kpc in order to be behind the warm HI at 
the same velocity. We analyze this quantitatively below. 



3.2 Is the amount of cold HI gas in the molecular 
cloud enough to produce significant HISA or 
HINSA against warm HI background? 

First we consider the conditions for HISA. The intensity 
of HISA is given by AT HISA = (T B (v) - T s )(l- e - T ). With 
a normal background HI temperature of about 100k, to 
produce a measurable H ISA signal at about 10 K level 
i|Anderson fc Banial 2009) requires an optical depth of HISA 
above 0.1. Such optical depth is produced by a cold HI col- 
umn of lxlO 20 cm -2 according to N H/ = 1.82x W 19 T s r Av, 
assuming the cold gas has a velocity dispersion of 5 km 
s _1 . The continuum brightness temperature of the south- 
ern bright spot of CTB109 is 30K, so the absorption feature 
at —51 km s _1 should be produced by cold HI clouds at 
—51 km s" 1 (with T s ^ 30K). From the absorption spectra, 
the column density of the HI absorbing clouds is 3x 
cm" 2 , i.e. 3xl0 20 cm -2 when T 3 = 10 K. This is more than 
enough to produce measurable HISA if the cloud is in front 
of the warm HI emission (100 K). 

Next we calculate the optical depth of an HI absorption 
line which is required to produce a HINSA feature with peak 



intensity of 3 a. The rms of the HI emission data is about 1.7 
K per channel calculated fr o m the HI spectrum of G108.77- 
0.95 (Fig. 2). iTavlor et ail l|2003l ) give that the rms noise 
level in brightness temperature for an empty channel ranges 
from 2.1 K to 3.2 K throughout the survey. We take the 
value of 2.1 K here. Based on a three component radiative 
transfer calculation and the assumption tha t the HI emission 
is opt ically thin, we have (Equation 8 of iLi fc Goldsmithl 
(|2003t ). 



T ab = [pT H i + T c - T x ](l - e 



(1) 



where the depth of absorption T a b is determined by the frac- 
tion of HI in the background p, the emission temperature 
Thi, the continuum background T c , and the cold HI ex- 
citation temperature T x . Assuming a uniform Galactic HI 
disk with a 17 kpc radius, it is estimated that p = 0.815 
toward the direction of CTB 109. Using the reasonable as- 
sumption of T c = 3.5K (including CMB and interstellar 
radiation field see ILi fc Goldsmithl l|2003h ). T H i = WOK 
and T x = 10K, an optically depth of r = 0.084 is required 
to produce an absorption line of 6.3 K peak depth (3 a). 
For a FWHM line- width of 2.0 km s" 1 of HINSA, this cor- 
responds to an upper limit of the column density of cold 
HI N H i = 3.2xl0 18 cm' 2 . Theoretically, the abundance of 
atomic H I inside molecular clou ds is in the range of 0.05 
to 0.27% l|Goldsmith fc Lill2005l ). Observ ationally, using an 
improved technique to measure the ratio, iKrco et all ( 2008) 
obtained a precise value for two clouds, i.e. 0.09% for L134 
and 0.53% for L1757. Taking the theoretical mean HI/H2 ra- 
tio of 0.16%, the non-detection of HINSA toward the large 
molecular cloud interacting with CTB 109 puts a 3<r upper 
limit on the cloud H2 column density at 2.0 x 10 21 cm -2 . 

The 12 CO spectrum from box 4 (Fig. 3) gives Wi2 CO = 
12 Kkms' 1 for the cloud. The 12 CO to H 2 conversion factor 
X is in the range of 2.2 to 2.8 x 10 20 [cm~ 2 /(Kkms' 1 )} base d 
on three calibration techniques (ISolomon fc Barrett! Il99ll ). 
This gives an H2 column density of 3 xlO 21 cm -2 for this 
cloud (X = 2.5xl0 20 cm" 2 /(Kkms' 1 )). This is enough to 
produce measurable HINSA at 4.5a" level if the cloud is in 
front of the warm HI at the same velocity. 

Fig. 3 (third row) shows a good example that a HISA 
feature at —50 km s" 1 is clearly detected nearby CTB 109. 
The CO cloud at —50 km s" 1 has low brightness tempera- 
ture of 0.3 K and shows no physical relation with the rem- 
nant. The Wi2 CO from the cloud is about 1 K km/s so it 
has an H2 column density of 2.5 xl0 20 cm -2 . This is about 
30 times lower than that from the cloud complex associated 
with Sh 152. This indicates that the cold HI cloud at —50 
km s" 1 is likely responsible for the absorption feature. 



3.3 Distance of 4 kpc to CTB 109/AXP IE 
2259+586 

The non-detection of both HISA and HINSA leads to the 
conclusion that the cloud at —55 km s" 1 is behind the warm 
HI and thus CTB 109, interacting with the cloud, is behind 
the velocity reversal in the Perseus arm. As noted above, the 
HI m odel distance of 4 kpc has a uncertainty of 20% (i.e. 0.8 
kpc). lBrand fc Blitj (|l993l) have studied the observed veloc- 
ity field of the outer galaxy and found residuals from circular 
motions. They found that these velocity residuals are con- 
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sistent with streams motions due to spiral arms, supporting 
the Roberts (and thus FM06) models. 

K02 estimated a distance for CTB 109 of 3±0.5 kpc 
based on two arguments. First, the HI absorption at —45 
km s" 1 in the extragalactic source is not seen in the spectra 
of CTB 109 or Sh 152, so they argued CTB 109 and Sh 152 
must be at the nearer distance in the spiral shock region. 
However, the extragalactic source is not near to CTB 109 
(outside of the area shown in Fig.l) and thus spatial varia- 
tions in the HI distribution can result in lack of significant 
HI at -45 km s" 1 in front of CTB 109 and invalidate that 
conclusion. Secondly they considered the spectroscopic dis- 
tances and radial velocities of 11 HII regions which in the 
Perseus arm nearby on the sky to CTB 109. These HII re- 
gions have an average distances of 3 kpc. However the two 
HII regions, Sh 152 ([108.77,-0.95]) and Sh 149 ([108.34, - 
1.12]), closest to CTB 109 on the sky have spectroscopic 
distances of 3.6±1.1 kpc and 5.4±1.7 kpc, respectively. This 
does not favour the small distance to CTB 109 of 3 kpc. 
These two HII regions have respective radial ve locities of 
-50.4 +0.5 km s" 1 (or -49.1±0.9 k m s" 1 from iLockmanl 
ljl989h and -53.1±1.3 km s" 1 from iBrand fc Blitj (|l993l ) 
which are similar to the radial velocity ( — 55 km s" 1 ) of the 
molecular cloud complex with which CTB 109 is interact- 
ing. This also support the larger distance for CTB 109 of 
4.0 kpc. 

CTB 109 is host to AXP IE 2259+586. DvK06 excluded 
small distance of 3 kpc for IE 2259+586 due to absence of 
red clump stars sufficiently highly reddened to be consistent 
with IE 2259+589's column density. There are two jumps (at 
Ay 3 and 8) in their Fig. 7 which shows a redding-distance 
relation along the line-of-sight to AXP IE 2259+586. They 
favour a large distance of 6 kpc for the AXP based on the 
large extinction of 8. The HI absorption spectrum of CTB 
109 (Fig. 2) excludes a large distance. Eg for 6 kpc the ra- 
dial velocity of FM06 model is about —80 km s , but no 
absorption is seen beyond —56km s" 1 . Our distance mea- 
surement of 4.0 kpc to AXP IE 2259+586 puts the AXP 
at the far edge of the Perseus arm and would give a small 
extinction of Ay ~ 3 for the AXP. The reddening stays at 
Ay ~ 3 out to about 6 kpc (Fig. 7 of DvK06), implying 
excess local extinction for the AXP at 4 kpc. Using Nj? = 
Ay x 2.3xl0 21 cm -2 , with inferred excess AAy = 5, gives 
excess ANh ~ 10 22 cm" 2 . This is quite plausible given that 
AXPs should be produced in SN explosions of massive stars 
and these SN also are known to produce large amount of 
dust in their ejecta. 

The new distance of 4 kpc to the AXP leads to a larger 
explosion energy of L 3X10 ergs than that based on 3 kpc 
l|Vink fc Kuiperll2006h . However this is not a large enough 
increase to argue that CTB 109 had an unusually large 
explosion energy. The same si tuation occurs for SN R Kes 
73/AXP IE 1841-045 system dTian fc Leahy||2008bl ). Both 
cases have larger distances than previously estimated but 
still give n ormal explosion energi es (10 51 ergs). This adds 
support to IVink fc Kuiperl ||2006 :)'s argument against the 
possibility that the magnetars are formed from rapidly ro- 
tating (less than a few millisecond) proto-neutron stars. We 
have applied a new technique for determining the distance to 
CTB 109, i.e. absence or presence of HISA and HINSA. This 
is a new useful tool for finding distances to other objects in 
the Galactic plane. 
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